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Abstract 

In this paper, we propose a novel approach for texture analysis based 
on artificial crawler model. Our method assumes that each agent can 
interact with the environment and each other. The evolution process con- 
verges to an equilibrium state according to the sot of rules. For each 
textured image, the feature vector is composed by signatures of the live 
agents curve at each time. Experimental results revealed that combining 
the minimum and maximum signatures into one increase the classification 
rate. In addition, we pioneer the use of autonomous agents for character- 
izing silk fibroin scaffolds. The results strongly suggest that our approach 
can be successfully employed for texture analysis. 
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1 Introduction 



The silk fibroin is a protein extracted from cocoons of Bombyx mori silkworm. It 
has been widely used in biomedical applications due to its high capacity to suture 
tissues, to regenerate bones and its biocompatibility with several types of human 
cells used in prosthesis PQ[2]. In the last years, researchers have proposed to 
improve the nanoscale features of silk fibroin by adding glycerol [2] ■ Though the 
presence of glycerol can provide better material properties, it also can alter the 
silk fibroin molecules interactions, damaging the result in its surface. Therefore, 
texture analysis methods emerge as a powerful tool for determining the suitable 
concentration of glycerol. 

Many methods for texture description have been proposed in the literature. 
These methods are based on statistical analysis of the spatial distribution (e.g., 
co-occurrence matrices [3j, local binary pattern j3] and entropy j5j), stochastic 
models (e.g., Markov random fields [6j), spectral analysis (e.g., Fourier descrip- 
tors [7], Gabor filters [H [9] and wavelets transform [10] ). complexity analy- 
sis (e.g., fractal dimension [llj). agent-based model (e.g., deterministic tourist 
walk [Tallin]). Although there are effective texture methods, they do not capture 
the richness of patterns of the silk fibroin scaffolds. 

In this paper, we present a methodology for classifying surface properties 
of silk fibroin by means of texture analysis. The texture description approach 
proposed here is based on the artificial crawler model |14U15j . We propose a new 
rule of movement that not only moves artificial crawler agents toward higher 
intensity, as well as to lower ones. We confirm that this strategy increases the 
discriminatory power and outperforms the state-of-the-art method. 

This paper is organized as follows. Section [2] details the original artificial 
crawler model. Section[3]presents our approach to characterize textured images. 
Section |4] discusses the results of the experiments. Finally, conclusions are given 
in Section [Sl 

2 The Original Artificial Crawler Model 

The first artificial crawler ( ACrawler) model was developed in [21 [13] . Let us 
consider that an image is a pair (T,/), consisting of a finite set T of pixels, 
and a mapping / that assigns to each pixel p — {xp,yp) S T an intensity I(p) 
ranging from to 255. A pixel of the intensity map / holds a neighborhood set 
ri{p) of pixels q, which d{p, q) < \f2 is the Euclidean distance between pixels p 
and q. Thereby, we assume the eight-connected neighbors. 

The original artificial crawler model assumes that each agent lives on one 
pixel of the image. At each time an agent i is characterized by two attributes: 
(1) e* holds a level of energy and (2) p\ occupies a spatial position in the im- 
age. First, n agents are born with identical energy e. Such energy can either 
wax or wane their lifespan according to energy consumption and infiuence of 
the environment. On images, the environment is treated as a 3D surface with 
different altitudes that correspond to grey values in z-axis of the images. Higher 
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intensities supply nutrients to the agents, while lower altitudes correspond to 
the land. Figure [l] shows a textured image and the peaks and valleys where the 
agents can live. 




(a) (b) 



Figure 1: The environment of the agents. On the most left is shown a textured 
image (a) and its respective 3D surface (b). 

The algorithm consists of a set of rules that comprises the evolution process: 

1. Born: Each agent i is born with the same energy e. 

2. Survival threshold: An agent i dies if its energy is below the threshold. 

Vt^i, if e* < Cmin then i dies 

3. Movement: 

{p*, if (a) is satisfied 
pLax^ if (b) is satisfied 
p^, if (c) is satisfied 

(a) Agents settle down if the grey level of its 8-neighbors are lower than 
itself. 

(6) Agents move to a specific pixel if there exist one of its 8-neighbors 
{Pmax) with higher intensity 

(c) If there exist more than one neighbor with higher intensity, an agent 
moves to the pixel that already was occupied (p^). 

4. Energy consumption: Each time t consumes one unity of energy. 
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5. Law of the jungle: An agent with higher energy eats up another with 
lower one. 

^ij ■ P^'t^ = P]'^\emax{el,e*} = niax{e*,e*} 

6. Gain of energy: It up to dates the energy absorption from the environ- 
ment, where A is a rate of absorption over the pixel I{p\)- 

V„e*+i = e* + A/(p*) 

7. Limit of energy: It bounds the maximum energy Bmax- 

V, ■•(=>*.> p = p 

*J ' j — ^maxi^j — ^max 

Agents that were born in lower altitudes areas can die in the evolution pro- 
cess, while individuals that reached to settle down in areas of higher altitudes 
have higher likelihood to remain alive. To quantify the multi-agent system, a 
curve of live agents at each time is obtained: 

^=[^(0),^(l),...,V(imax)] (1) 

where ■(/'(t) is the number of live agents at time t and tmax is maximum iteration. 



3 A Novel Approach with Artificial Crawler to 
Texture Analysis 

The artificial crawler model described above consists of moving agents to a 
neighbor pixel with the highest intensity. Despite the promising results, this 
idea does not extract all the richness of textural pattern. Our approach differs 
from the original ACrawler model in terms of movement: each agent is not only 
able to move to the higher altitudes as well as to lower ones. It allows the model 
to extract the details present in peaks and valleys of the images. 

First, the agents move to higher intensities as the original artificial crawler 
method. Thus, the artificial crawlers are performed using this rule and the 
curve ^max is obtained. Throughout the paper, this rule of movement will be 
referred as max. We can observe that the original artificial crawler method only 
models the peaks of a textured image. To obtain a robust and effective texture 
representation, we propose to move artificial crawlers toward lower intensities 
— this rule of movement will be referred throughout the paper as to min. In 
our approach, artificial crawlers are randomly placed in the image with initial 
energy e. Then, the movement step is modified as follows: 

Vi :e* >e„i„,p*+i = /(p*) 



p\, if (a) is satisfied 
f{p) = { Plnini if (b) is satisfied 
p\^, if (c) is satisfied 
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(a) Agents settle down if the grey level of its S-neighbors are higher than itself. 
(6) Agents move to a specific pixel if there exist one of its 8- neighbors [plnin) 
with lower intensity 

(c) If there exist more than one neighbor with lower intensity, an agent moves 
to the pixel that already was occupied {pin)- 

The multi-agent systems using the rule of movement min is characterized as 
the original method by using the number of live agents at each time. Considering 
that now we have two rules of movement, the final feature vector of our approach 
is composed by the concatenation of fmax and ipmin according to: 

<P = [y^max, (Pmrn] (2) 

Figure [2] shows the curves of the evolution process. We took two classes of 
textures (on the top right-hand corner in Figure [2| from the album the Brodatz 
[l6] to illustrate the separability. On the left. Figure [5] shows the number of 
live agents using the rule of movement min, while the curve for the rule of 
movement max is shown on the right of Figure [2] The experimental results 
below corroborate the importance of both rules of movement in the texture 
modeling. 




5 1 1 5 20 25 30 35 40 5 1 1 5 20 25 30 35 40 

Iteration Iteration 



Figure 2: Curve of live agents using the rules of movement (left) max and 
(right) min. 



4 Experimental Results 

In this section, we demonstrate the effectiveness of our approach. We first 
outline details of the experimental setup, and then, experiments carried out 
on two datasets are discussed: Brodatz and silk fibroin. Here is described the 
whole process for image acquisition of silk fibroin scaffolds. Besides, we show 
comparative results with different texture methods. 
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4.1 Experimental Setup 

The proposed method was first evaluated over texture classification experiments 
by using images extracted from Brodatz album |16| . This album is considered a 
well-known benchmark for evaluating texture recognition methods. Each class 
is composed by one image divided into nine sub-images non-overlapped. A total 
of 440 images grouped into 40 classes was considered. Each image has 200 x 200 
of size and 256 grey levels. One example of each class is shown in Figure |3] 




Figure 3: Example of 40 Brodatz texture classes used in the experiment. Each 
image has 200 x 200 pixels and 256 gray levels. 

The texture classification was carried out for ten-fold cross validation to 
avoid bias. At each round, we randomly divide the samples of each class into 
ten subsets of the same size, i.e., nine for training and the remaining for testing. 
The results are reported as the average value over the ten runs. For classification, 
we adopted the model Linear Discriminant Analysis (LDA). The underlying idea 
is to maximize the Euclidian distance between the means of the classes, while 
minimizing the within-class variance. For further information we refer to [T7]. 

Linear Discriminant Analysis (LDA) [18] was selected since it is well founded 
in statistical learning theory and has been successfully applied to various object 
detection tasks in computer vision. LDA, originally proposed by Fisher, com- 
putes a linear transformation (T S 3?''^") of D, which D € 3?'^^" is a matrix 
and d denotes the number of features and n number of samples. 

We optimized two parameters of the artificial crawler model: the number of 
agents and the way that agents move in the evolution process. The number of 
agents placed on the pixels of the image was inital set to 1000 with a coverage 
rate of 10%, varying from 1000 to 35000. In our experiments, all agent was born 
with an initial energy e of 10 units and the loss for each iteration consumes 1 
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unit of energy. The absorption rate was set to 0.01 in terms of the current pixel. 
For the survival threshold and the upper bound of energy were set to 1 and 12 
units, respectively. 



4.2 Performance Evaluation 



Experiment 1: First, we perform an analysis of our method on the Brodatz 
dataset. Figure [4]^a) presents the correct classification rate versus the number of 
iterations. The results for the original artificial crawler is shown as curve max 
while the results for our method is shown as curve min U max. For a complete 
comparison, we also provide the results for an method which agents move to 
pixels with lower intensity — curve min. As can be seen, the proposed method 
provided the highest correct classification rates for all values of iterations. These 
experimental results indicate that the proposed method significantly improves 
performance over the traditional methods. We can also observe that the rule 
min provided higher rates than the rule max, given the idea that valleys are 
more discriminative than peaks in the Brodatz dataset. 

Another important parameter of the artificial crawler methods is the number 
of agents. Figure|4][b) shows the correct classification rates versus the number of 
agents. As in the previous experiment, our method achieved the highest rates 
compared to the other two strategies. Again, the rule min provided higher 
rates than the rule max. Another important observation from Figure |4]jb) 
is that using a few agents, the methods achieved good classification results, 
which makes the artificial crawler methods suitable for real time applications. 
Using these two plots, we can determine the best parameters of our method to 
tmax — 41 and n — 27k. 




Figure 4: Comparison of artificial crawler methods for different values of (a) 
iterations and (b) number of agents in the Brodatz dataset. 

The results of the proposed method are compared with existing texture anal- 
ysis methods in Table [ij It is observed that the our method outperforms the 
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state-of-the-art. The highest classification rate of 98.25%(±1.69) was obtained 
by our method, which is followed by a classification rate of 95.25%(±3.43) ob- 
tained by the Gabor filter, one of the most traditional texture analysis method. 



Method 


Images correctly classif. 


Correct classification(%) 


Fourier descriptors [7] 


346 


86.50 (±6.58) 


Co-occurrence matrices [3] 


365 


91.25 (±2.65) 


Original artificial crawler [14] 


372 


93.00 (±5.50) 


Gabor filter 9 


381 


95.25 (±3.43) 


Proposed method 


393 


98.25 (±1.69) 



Table 1: Experimental results for texture methods in the Brodatz dataset. 



Experiment 2: In this experiment, we present a comparative study of our ap- 
proach to assess the quality of the silk fibroin scaffolds. Our goal is to provide an 
effective method to support the visual analysis, thus reducing the subjectiveness 
of conclusions based on the human analysis. The potential of the silk fibroin is 
enhanced by including glycerol solutions during scaffold formation [2]. In gen- 
eral, such concentration can range from 0% to 10% with step of 2.5%. As far as 
the authors know, this paper is the first to report a method for characterizing 
the silk fibroin scaffolds. This dataset contains 5 classes, each of 10 200 x 200 
pixel images. Figure [s] shows three samples for each concentration. 
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Figure 5: Samples for each glycerol concentration. The first column corre- 
sponds to 0% of concentration, the second 2.5%, and so on up to 10%. 

We perform the same experiment to determine the best parameters of iter- 
ation and number of agents in the Silk Fibroin dataset. Figure |6|^a) presents 
the evaluation of t while Figure [6]jb) presents the evaluation of n for different 
artificial crawler methods. Using both plots, we found that the best results are 
achieved for tmax — 7 and n — 28k. 

In the silk fibroin dataset, our method achieved highest classification rates 
when compared with traditional texture analysis methods. The experimental 
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Figure 6: Comparison of artificial crawler methods for different values of (a) 
iterations and (b) number of agents in the silk fibroin dataset. 



results, presented in Table [2] shows that our method achieved a classification 
rate of 96%(±8.43). These experimental results indicate that our method is 
consistent and can be applied in real-world applications. 



Method 


Images correctly classif. 


Correct classification(%) 


Fourier descriptors [7] 


39 


78.00 (±22.01) 


Co-occurrence matrices [3] 


47 


94.00 (±9.66) 


Original artificial crawler [14] 


42 


84.00 (±15.78) 


Gabor filter |9] 


31 


62.00 (±19.44) 


Proposed method 


48 


96.00 (±8.43) 



Table 2: Experimental results for texture methods in the silk fibroin dataset. 



5 Conclusion 

In this paper we presented a novel approach based on artificial crawler for tex- 
ture classification. We have demonstrated how the feature vector can be im- 
proved by combining min and max curves, instead of using only the strategy for 
the maximum of intensity of the pixels. Although traditional methods of tex- 
ture analysis have provided satisfactory results, the approach proposed here has 
proved to be superior for characterizing textures. We have tested our proposal 
on the most popular benchmark for texture analysis and find that it produces 
good classification results. Furthermore, we tested our approach on the silk 
fibroin scaffolds analysis and results indicate that our method is consistent and 
can be applied in real-world applications. 

The results support the idea that our approach can be used as a feasible step 
for many analysis not only applications on tissue engineering. In addition, our 
results can be improved by studying the setting of our method. We have already 
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implemented other variation. As part of the future work, instead using random 
sampUng of agents we plan focus on evaluating the deterministic sampling, i.e, 
each pixel of the image is initialized with on agent. 
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